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1 Introduction 

The abundance of dark matter haloes is now believed to be a sensitive cosmological test ( [1] ). 
The halo mass function (hereafter MF) is used as an ingredient for the clustering models of 
haloes and galaxies ([2], [3]). Also it is the subject of the “satellite abundance” problem ([4]). 

The most popular, useful and developed analytical model of MF prediction is the Press- 
Schechter method ([5], [6], [7], [8]). It is based on the spherical collapse of dusty matter cloud 
which make this model one of the simplest. However, numerical simulations of Large Scale 
Structure show that the Press-Schechter model predicts somewhat wrong distributions ([9], 
[10], [11], [12], [13], [14]). There exists a more precise analytical method of modeling of the 
MF, the Sheth-Mo-Tormen ([15]) method, in which a spherical collapse is substituted by an 
ellipsoidal collapse. This gives rise to more accurate results especially at high masses. The 
results of simulated MFs are usually represented as analytical hts ([16], [17]) which nowadays 
have the precision of 5-10%. There also exist more sophisticated methods of MF calculation 
from the initial random density field without running a simulation, e.g. the PINOCCHIO 
code ([18]). 

The most reliable results on the MF are now obtained with the help of large N-body 
simulations. However they are costly, especially when the goal is to investigate the depen¬ 
dence of the MF on cosmological parameters. On the other hand, nowadays, in the era of 
precise cosmology, the range of parameter change is quite small, e.g. the difference in the 
matter density Qm between WMAP5 and Planck cosmologies is less then 0.05. 

It is therefore attracting to scale the results of existing N-body simulations instead of 
doing new ones. In order to do this one needs to calculate analytically how the MF changes 
with the small change of cosmological parameters. This calculation can be done basing on 
the simplest Press-Schechter method. Examples of such scaling can be found in [19], [20], [21] 
who used a simulated MF for Gaussian initial conditions and multiplied it by a correction 
factor extracted from the ratio of a non-Gaussian Press-Schechter MF to a Gaussian one. A 
similar approach was also used by [22] and [23] who computed a MF in different environments 
by taking a high precision global MF (Sheth-Mo-Tormen in [22] and a simulated one in [23]) 
and multiplying it by a bias factor computed from an excursion set approach. 

Despite the fact that the rescaling of MF has been used in several papers, there is no 
exhaustive investigation of the accuracy of this approach in literature up to date. E.g. [23] 
computed the number of haloes in the mass range 10^ — 10 ®Mq as a function of large-scale 
density and found that using the rescaling of either N-body MF or Sheth-Mo-Tormen MF for 
this purpose is much more accurate than using linear bias model with Press-Schechter MF. 
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Figure 1. Halo mass functions for two sets of cosmological parameters: flm = 0.2 and flm = 0.3 
obtained with Press-Schechter method and simulations. 


The rescaling of MF may turn out to be useful for the exploration of a variety of non¬ 
standard cosmologies as well as fine-tuning the cosmological parameters to fit the observed 
Universe. This motivated us to investigate the precision of the rescaled mass function in this 
paper. For this purpose we run a set of simulations with varying cosmological parameters: 
matter density ^Im, amplitude of the power spectrum as, slope of the initial power spectrum 
Us- The accuracy of the rescaled MF is then expressed as a function of the difference of these 
parameters between the source and target cosmologies. 

There is another way of rescaling: treating the MF as a function of a{M), the RMS of 
matter density fluctuations on scale M. Press-Schechter theory implies that the dependency 
of MF on a[M) does not depend itself on the cosmological parameters, i.e. it is universal. 
During the past decade this universality has been intensively checked ([11], [16], [17], [24]). 
The conclusion is that the MF may differ from the universal one by 5-10% in case of Friends- 
of-Friends haloes, or 20-50% for spherical haloes ([17]). 

Another technique of rescaling was demonstrated by [25]. They created the algorithm 
that allows to scale the output of a cosmological N-body simulation carried out for one specific 
set of cosmological parameters so that it faithfully represents the growth of structure in a 
different cosmology. 

The paper is organized in the following way. Firstly we briefly outline the theoretical 
background on which the rescaling method is based. In section 2 we present results and test 
the accuracy of the method for different cosmological models. In the last section we present 
summary of our work and discuss some astrophysical implications. 

2 Method description and tests 
2.1 Method of MF rescaling 

Firstly, let us describe the basics of the method with a simple example. In figure 1 we present 
two mass functions produced in simulations and two mass functions calculated using Press- 
Schechter method for two different sets of parameters: in one case Dm = 0.2 (Da = 0.8, 
hereafter we will omit Da implying flat cosmology) and in another one Dm = 0.3. 
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Figure 2. The ratio of two Press-Schechter mass functions and two simulated mass functions. 


The values flm = 0.2 and VLrn = 0.3 were chosen for demonstration purposes. It is 
seen that the curve shapes are similar for theory and simulations (points near the right edge 
contain large errors because of small number of haloes in this region). The similarity is 
demonstrated even more clearly in the figure 2, where the ratio of two Press-Schechter MFs 
and two simulation MFs are shown. 

Due to this finding we suppose that the ratio between pairs of simulated and theoretical 
MFs at any given mass is equivalent and corresponds to the next equations (the dependency 
of MFs on mass is omitted): 

The ratio between Press-Schechter MFs: 


PS{n^ = 0 . 2 ) 

PF(D^ = 0.3)' 


The ratio between simulated MFs: 


_ F/M(D™ = 0.2) 
SIM{n^ = 0.3) ■ 


From the similarity of the ratios: 


Rps ~ Rsim, 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


it follows that the new MF can be expressed as: 

Resale,^ 0 . 3 ) ^ ^ 

In figure 3 we present an original simulated mass function for = 0.3 and an extrap¬ 
olation of the mass function for this set of parameters from Dm = 0.2 using (2.4). They are 
extremely close. The error is very small mostly everywhere except for high masses. The main 
goal of our research was to quantify this error as a function of the difference in cosmological 
parameters. 
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Figure 3. The MF from the original simulation for ftm = 0.3 and the reconstructed MF by using 
(2.4) for rtm = 0.3. 


2.2 Simulations 

To test the accuracy of rescaling we run a set of N-body simulations using GADGET-2 code 
([26]) with dark matter only. Volume of all the simulations is a cubic box of 500/i“^ Mpc size, 
the total number of particles is 512^, the gravitational softening parameter is 25h~^ kpc. For 
the initial conditions we used the N-GenIG code by V. Springel. We held the phases of the 
initial velocity field the same for all the simulations. The cosmological model is flat Universe 
with ACDM -cosmology where = 1- The simulations were run on the supercomputer 

of the Nuclear Physics and Astrophysics Division of P.N. Lebedev Physical Institute. 

For finding the haloes we used the Amiga’s Halo Finder (AHF) ([27], [28]). The mini¬ 
mum number of particles in a halo is set to 60, the total number of haloes in each simulation 
is about 10^ at z = 0. For the MF construction the Friends-of-Friends ([29]) method is used 
quite often e.g. in [16], so we also used RockStar halo finder ([30]). The virial overdensity A 
is calculated for each set of cosmological parameters using spherical top-hat collapse model 
([31]). Only the isolated haloes are used. 

We vary three cosmological parameters, Dm, <78 and n^. Our set of simulations is sum¬ 
marized in Table 2.3. In addition to our own simulations we use two public available simula¬ 
tions: Bolshoi and BolshoiP ([32]). The halo catalogues were obtained from the COSMOSIM 
database (www.cosmosim.org). The haloes were identified using the standard overdensity 
criterion with A = 360 • Pback with the spherical overdensity-based BDM halo hnder. 

2.3 The testing of the method accuracy 

In order to measure the accuracy of the method we use the maximum difference between two 
cumulative MFs. In simulations and Press-Schechter method we obtain differential MFs (here 
the division into 1000 bins is used), integrate them and then interpolate for the same mass 
interval. Finally, we take a maximum of the difference between the original and calculated 
MFs divided by the original MF. The maximal value defined this way we call the error of the 
method and it is shown in figures 4-6. For each graph there are presented the errors for the 
rescaled MF, for pure Press-Schechter, Sheth-Mo-Tormen and Tinker fits. 
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Name 


Ha 

0-8 

Us 

Sim(nm=0.2) 

0.2 

0.8 

0.8 

1.0 

Sim(nm=0.25) 

0.25 

0.75 

0.8 

1.0 

Sim(I2m=0.3) 

0.3 

0.7 

0.8 

1.0 

Sim(I2m=0.35) 

0.35 

0.65 

0.8 

1.0 

Sim(nm=0.4) 

0.4 

0.6 

0.8 

1.0 

Sim((T8=0.9) 

0.25 

0.75 

0.9 

1.0 

Sim((T8=1.0) 

0.25 

0.75 

1.0 

1.0 

Sim(ns=0.9) 

0.25 

0.75 

0.8 

0.9 

Sim(ns=0.8) 

0.25 

0.75 

0.8 

0.8 

Bolshoi 

0.27 

0.73 

0.82 

0.95 

BolshoiP 

0.30711 

0.69289 

0.82 

0.96 


Table 1. Parameters of simulations used in tests 
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Figure 4. The accuracy of MF calculation methods in dependence upon the Qm difference. 


From the figures 4-6 we see that for the small change in cosmological parameters the 
method of rescaling has much better accuracy than Press-Schechter, Sheth-Mo-Tormen and 
Tinker MFs. The accuracy of the rescaled MF also can be quantihed in the following way: 


AMF = 36% AQrn, 

(2.5) 

AMF = 18% An,, 

(2.6) 

AMF = 30% Afj8. 

(2.7) 


The result of test of the rescaling method with the Bolshoi and BolshoiP simulations 
is shown in figure 7. In it the original MF for the Planck parameters together with the MF 
obtained using equation (2.4) from WMAP-based simulation are shown on the top. On the 
bottom the ratio of MFs is shown. The net error defined above is 6% for these simulations. 
It is 4 times higher than expected from (2.5), however still much less than for pure Press- 
Schechter (70%) or Sheth-Mo-Tormen (20%). One should note that Bolshoi and BolshoiP 
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Figure 5. The accuracy of MF calculation methods in dependence upon the as difference. 
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Figure 6. The accuracy of MF calculation methods in dependence upon the Us difference. 


cover two orders of magnitude larger mass range than our own simulations, so the error in 
the cumulative MF is being accumulated. If we consider only haloes with M > in 

Bolshoi and BolshoiP, the error decreases to 2.8%. 

From figure 7 it is seen that at lower masses, < 3 • the rescaled MF, Sheth- 

Mo-Tormen and Tinker fits all give the same error: they underpredict the number of haloes by 
several percent. At larger masses Sheth-Mo-Tormen and Tinker fits overpredict the number 
of haloes while the rescaled MF is the most precise. 

We also test how our results depend on the halo finder used (FoF or spherical over¬ 
density) and find the impact of halo identification algorithm is quite small: the error of the 
rescaled MF is about 20% larger for the FoF haloes than for the spherical overdensity. This 
was also noticed by [33] who found that the halo bias factor is insensitive on how the haloes 
are identified. 
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Figure 7. Top: the original MF for the BolshoiP simulation and the MF calculated by rescaling 
using equation (2.4) and the MF for the Bolshoi simulation. 

Bottom: the ratio of different MF approximations to the MF measured from BolshoiP simulation: 
rescaled MF (solid line with error bars), Press-Schechter (dashed line), Sheth-Mo-Tormen (dash-dot 
line), Tinker fit (dash-dot-dot). 

3 Summary and discussion 

We have investigated the accuracy of the method of calculating a MF of dark matter haloes 
based on the rescaling of a MF measured in simulations. The method is described by equation 
(2.4). Our tests show that when the difference of cosmological parameters between the 
target and the source cosmologies is small, the accuracy is very high, reaching few percent 
(see figures 4-7 and equations (2.5)-(2.7)). This method can be applied to problems which 
require varying cosmological parameters, such as fitting of observational data. 

In comparison to the most accurate existing fits to simulations (e.g. [17]), the rescaled 
MF has several advantages. Firstly, it predicts the number of haloes for the same technique 
of halo identification which is used in the source simulation, which may differ from one used 
in the available fits. Secondly, we expect that it can be used for non-standard cosmological 
models, e.g. with warm dark matter, non-standard dark energy, etc. 

Another potential application of the rescaling may be the increase of precision on the 
high mass end. The idea is to run an existing simulation past z = 0 in order to facilitate 
the growth of high mass haloes and thus increase statistics. Then the MF can be rescaled 
back to z = 0 using equation (2.4). For example, in order to reduce the statistical errors two 
times one needs to increase the number of haloes four times. For masses M > 10^^Mq to do 
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so we need to go in future by z = —0.33 or 6 Gyr (in PLANCK cosmology). 

One can also try to use the universal MF approach for a similar procedure of rescaling, 
i.e. to express a MF from a source simulation as a function of a and then to substitute a{M) 
for the target cosmology. In this approach the mass range of the resulting MF differs from 
that of the source simulation, what will not allow to use this method to increase the precision 
at high mass end as described above. Another problem will arise in the case of WDM since a 
weakly depends on M for masses smaller than the filtering mass. Therefore M[a) becomes 
undefined for large enough a and, hence, this approach may fail in some cases when the 
source cosmology is WDM. 

Acknowledgment s 

The authors thank Mikhail Zavertyaev for providing computational time and tuning the 
NPAD LPI cluster. We also thank Kristin Riebe for helping with the Cosmosim database, 
Ravi Sheth and Cristiano Porciani for helpful discussion. The work of S. Pilipenko is sup¬ 
ported by Basic Research Program P-7 of the Presidium of the Russian Academy of Sciences, 
by “UNK FIAN” and by grant of President of Russian Federation for Support of Leading 
Scientifics School NSh-4235.2014.2. VY was partly supported for this research by the Dy¬ 
nasty Foundation and partly supported through a research contract from the International 
Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities 
of Bonn and Cologne. 

The CosmoSim database used in this paper is a service by the Leibniz-Institute for 
Astrophysics Potsdam (AIP). The MultiDark database was developed in cooperation with 
the Spanish MultiDark Consolider Project CSD2009-00064. The Bolshoi and MultiDark 
simulations have been performed within the Bolshoi project of the University of California 
High-Performance AstroComputing Center (UC-HiPACC) and were run at the NASA Ames 
Research Center. The MultiDark-Planck (MDPL) and the BigMD simulation suite have been 
performed in the Supermuc supercomputer at LRZ using time granted by PRACE. 

We also used HMFcalc (http://hmf.icrar.org/) by [34] to verify our calculations of Press- 
Schechter, Sheth-Mo-Tormen and Tinker fits. 

References 

[1] A. Vikhlinin, A. V. Kravtsov, R. A. Burenin, H. Ebeling, W. R. Forman, A. Hornstrup et ah, 
Chandra Cluster Cosmology Project III: Cosmological Parameter Constraints^ 

Astrophys. J 692 (Feb., 2009) 1060-1074, [0812.2720]. 

[2] A. Cooray and R. Sheth, Halo models of large scale structure, 

Physics Reports 372 (Dec., 2002) 1-129, [astro-ph/0206508]. 

[3] M. P. Viero, L. Moncelsi, F. Mentuch and et al.. Measuring star formation in high-z massive 
galaxies: a mid-infrared to submillimetre study of the GOODS NICMOS Survey sample, 

Mon. Not. R. Astron. Soc 421 (Apr., 2012) 2161-2169, [1008.4359]. 

[4] M. Boylan-Kolchin, J. S. Bullock and M. Kaplinghat, The Milky Way’s bright satellites as an 
apparent failure of KCDM, Mon. Not. R. Astron. Soc 422 (May, 2012) 1203-1218, 

[1111.2048]. 

[5] W. H. Press and P. Schechter, Formation of Galaxies and Clusters of Galaxies by Self-Similar 
Gravitational Condensation, Astrophys. J 187 (Feb., 1974) 425-438. 


- 8 - 


[6] J. R. Bond, S. Cole, G. Efstathiou and N. Kaiser, Excursion set mass functions for hierarchical 
Gaussian fluctuations, Astrophys. J 379 (Oct., 1991) 440-460. 

[7] R. G. Bower, The evolution of groups of galaxies in the Press-Schechter formalism, 

Mon. Not. R. Astron. Soc 248 (Jan., 1991) 332-352. 

[8] C. Lacey and S. Cole, Merger rates in hierarchical models of galaxy formation, 

Mon. Not. R. Astron. Soc 262 (June, 1993) 627-649. 

[9] C. Lacey and S. Cole, Merger Rates in Hierarchical Models of Galaxy Formation - Part Two - 
Gomparison with N-Body Simulations, Mon. Not. R. Astron. Soc 271 (Dec., 1994) 676, 
[astro-ph/9402069] . 

[10] R. K. Sheth and G. Tormen, Large-scale bias and the peak background split, 

Mon. Not. R. Astron. Soc 308 (Sept., 1999) 119-126, [astro-ph/9901122]. 

[11] A. Jenkins, C. S. Frenk, S. D. M. White, J. M. Colberg, S. Cole, A. E. Evrard et ah, The mass 
funetion of dark matter haloes, Mon. Not. R. Astron. Soc 321 (Feb., 2001) 372-384, 

[astro-ph/0005260] . 

[12] G. Efstathiou, C. S. Frenk, S. D. M. White and M. Davis, Gravitational clustering from 
scale-free initial conditions, Mon. Not. R. Astron. Soc 235 (Dec., 1988) 715-748. 

[13] G. Efstathiou and M. J. Rees, High-redshift quasars in the Gold Dark Matter cosmogony, 

Mon. Not. R. Astron. Soc 230 (Feb., 1988) 5P. 

[14] S. D. M. White, G. Efstathiou and C. S. Frenk, The amplitude of mass fluetuations in the 
universe, Mon. Not. R. Astron. Soc 262 (June, 1993) 1023-1028. 

[15] R. K. Sheth, H. J. Mo and G. Tormen, Ellipsoidal collapse and an improved model for the 
number and spatial distribution of dark matter haloes, 

Mon. Not. R. Astron. Soc 323 (May, 2001) 1-12, [astro-ph/9907024]. 

[16] M. S. Warren, K. Abazajian, D. E. Holz and L. Teodoro, Precision Determination of the Mass 
Function of Dark Matter Halos, Astrophys. J 646 (Aug., 2006) 881-885, [astro-ph/0506395]. 

[17] J. Tinker, A. V. Kravtsov, A. Klypin, K. Abazajian, M. Warren, G. Yepes et ah. Toward a 
Halo Mass Function for Precision Gosmology: The Limits of Universality, 

Astrophys. J 688 (Dec., 2008) 709-728, [0803.2706]. 

[18] P. Monaco, E. Sefusatti, S. Borgani, M. Crocce, P. Fosalba, R. K. Sheth et ah. An accurate tool 
for the fast generation of dark matter halo catalogues, 

Mon. Not. R. Astron. Soc 433 (Aug., 2013) 2389-2402, [1305.1505]. 

[19] M. LoVerde, A. Miller, S. Shandera and L. Verde, Effects of scale-dependent non-Gaussianity 
on cosmological structures. Journal of Gosmology and Astroparticle Physics 4 (Apr., 2008) 14, 
[0711.4126]. 

[20] X. Kang, P. Norberg and J. Silk, Can a large-scale structure probe cosmic microwave 
background-constrained non- Gaussianity ?, 

Mon. Not. R. Astron. Soc 376 (Mar., 2007) 343-347, [astro-ph/0701131]. 

[21] M. Grossi, K. Dolag, E. Branchini, S. Matarrese and L. Moscardini, Evolution of massive 
haloes in non-Gaussian scenarios, Mon. Not. R. Astron. Soc 382 (Dec., 2007) 1261-1267, 
[0707.2516]. 

[22] R. Barkana and A. Loeb, Unusually Large Fluctuations in the Statistics of Galaxy Formation 
at High Redshift, Astrophys. J 609 (July, 2004) 474-481, [astro-ph/0310338]. 

[23] K. Ahn, 1. T. lliev, P. R. Shapiro and C. Srisawat, Non-linear bias of cosmological halo 
formation in the early universe, Mon. Not. R. Astron. Soc 450 (June, 2015) 1486-1502, 
[1407.2637]. 


- 9 - 


[24] S. Bhattacharya, K. Heitmann, M. White, Z. Lukic, C. Wagner and S. Habib, Mass Function 
Predictions Beyond ACDM^ Astrophys. J 732 (May, 2011) 122, [1005.2239]. 

[25] R. E. Angulo and S. D. M. White, One simulation to fit them all - changing the background 
parameters of a cosmological N-body simulation, 

Mon. Not. R. Astron. Soc 405 (June, 2010) 143-154, [0912.4277]. 

[26] V. Springel, The cosmological simulation code GADGET-2, 

Mon. Not. R. Astron. Soc 364 (Dec., 2005) 1105-1134, [astro-ph/0505010]. 

[27] S. P. D. Gill, A. Knebe and B. K. Gibson, The evolution of substructure -LA new 
identification method, Mon. Not. R. Astron. Soc 351 (June, 2004) 399-409, 

[astro-ph/0404258] . 

[28] S. R. Knollmann and A. Knebe, AHF: Amiga’s Halo Finder, 

Astrophys. J., Supplement 182 (June, 2009) 608-624, [0904.3662]. 

[29] M. Davis, G. Efstathiou, C. S. Erenk and S. D. M. White, The evolution of large-scale structure 
in a universe dominated by cold dark matter, Astrophys. J 292 (May, 1985) 371-394. 

[30] P. S. Behroozi, R. H. Wechsler and H.-Y. Wu, The ROCKSTAR Phase-space Temporal Halo 
Finder and the Velocity Offsets of Gluster Gores, Astrophys. J 762 (Jan., 2013) 109, 
[1110.4372]. 

[31] M. A. K. Gross, Gdm-Variant Gosmological Models: Nonlinear Gomparisons Including 
Large-Scale Constrained Realizations. PhD thesis, UNIVERSITY OF GALIFORNIA, SANTA 
CRUZ, 1997. 

[32] A. A. Klypin, S. Trujillo-Gomez and J. Primack, Dark Matter Halos in the Standard 
Cosmological Model: Results from the Bolshoi Simulation, Astrophys. J 740 (Oct., 2011) 102, 
[1002.3660]. 

[33] H. J. Mo and S. D. M. White, An analytic model for the spatial clustering of dark matter 
haloes, Mon. Not. R. Astron. Soc 282 (Sept., 1996) 347-361, [astro-ph/9512127]. 

[34] S. G. Murray, C. Power and A. S. G. Robotham, HMFcalc: An online tool for calculating dark 
matter halo mass functions, Astronomy and Computing 3 (Nov., 2013) 23-34, [1306.6721]. 


- 10 


